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Abstract 

The stability of a thermocapillary flow in an extended cylindrical geometry is analyzed. This flow 
occurs in a thin liquid layer with a disk shape when a radial temperature gradient is applied along the 
horizontal free surface. Besides the aspect ratio, a second parameter related to the local curvature is 
introduced to describe completely the geometrical effects. We recover classical hydrothermal waves 
as predicted by Smith and Davis pQ but the properties of these waves are shown to evolve with the 
curvature parameter, thus leading to a non uniform pattern over the cell. Moreover, it is shown 
that the problem is not invariant with respect to the exchange of the hot and cold sides. 

Keywords: Hydrodynamic instabilities, hydrothermal waves, thermocapillary flow 

Resume 

Nous etudions la stabilite lineaire d'un ecoulement thermocapillaire en geometrie cylindrique e- 
tendue. Un tel ecoulement est produit dans un disque de fluide dont la surface libre horizontalc 
est soumise a un gradient de temperature purcmcnt radial. Outre lc rapport d'aspect, un sec- 
ond parametre lie a la courbure locale est introduit pour caracteriser la geometrie du probleme. 
L'instabilite en ondes hydrothermales predite par pQ est retrouvee mais les proprietes des ondes 
sont alterees par la courbure locale ce qui cxplique l'existence de structures non uniformes. La 
dissymetrie du probleme vis-a-vis d'une inversion du gradient de temperature est aussi mise en 
evidence. 

Mots cles: Instabilites hydrodynamiques, ondes hydrothermales, ecoulement thermocapillaire 
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1 Sketch of the problem 



When a horizontal temperature gradient is imposed over a thin fluid layer with a free surface, a basic 
flow is present due to the Marangoni effect, i.e., surface tension dependence on the temperature. 
When the temperature difference is small, the characteristic velocity of the flow is proportional to 
the thermal constraint and the flow is stationary. When the temperature difference is increased 
above a threshold value, the basic flow is unstable to propagating waves called hydrothermal waves 
and first predicted by Smith and Davis pp for a zero gravity environment. 

The stability analysis of the basic flow has been addressed by many authors in Cartesian co- 
ordinates 0, corresponding to experiments in rectangular cavities |3J, but a very few theoretical 
studies were devoted to the cylindrical geometry 0] though experiments are often conducted in this 
configuration [5] . The present paper aims to fill the gap and investigates the effect of curvature on 
the critical properties of hydrothermal waves. 

The section of a typical cell in a vertical diametrical plane is depicted on Fig. ^ We note h 
the depth of the fluid layer, and R\ , i?2 = Ri + L the radii of the internal and external cylindrical 
boundaries that confine the fluid in the horizontal direction. The horizontal temperature gradient 
is supposed to be applied as a temperature difference between the inner and the outer boundaries, 
considered as isothermal with temperature T\ and Ti respectively. The bottom of the cell is 
considered as perfectly conducting [o| and a mixed thermal boundary condition is used for the 
top free surface 2 . The fluid is assumed to be incompressible with kinematic viscosity thermal 
diffusivity k, density p, thermal expansion coefficient a. The surface tension a of the interface varies 
with the temperature T of the liquid and its derivative da / dT = —7, is considered as constant and 
negative (7 > 0). 



2 Formulation of the basic flow 

To fully characterize the geometrical configuration, two non-dimensional parameters are required. 
First, the aspect ratio e = h/L which is supposed to be small in the present extended geometry 
(s < 0.03). Second, the local curvature proportional to the inverse of the radius will be represented 
in non dimensional form by L/r. Its maximum value: cp — L/Rx, is taken as the additional 
parameter characteristic of the cylindrical geometry as a whole. When the change of variable: 
r = R\ + LX, is made the local curvature is expressed as 

*(*) = ^— = - 
V 1 \ + 4>X r 

The basic flow is governed by the Navier-Stokes equations in the Boussinesq approximation together 
with the heat equation. Owing to rotational invariance of the system about the vertical axis e z of 
the cell, the basic flow has no azimuthal velocity component. The governing equations are written 
in non-dimensional form by taking the following scales : h for the vertical coordinate z, v/h for 
vertical velocity w and radial velocity u, h2/v for the time t, vIp^LfKi for the pressure P and 
T 2 — Ti for the temperature difference T — T\ . 

The boundary conditions represent a conducting rigid bottom at z — and a mixed thermal 
condition on the top free surface at z = 1 : 

u = w = 0, and T — T C (X) on z = 0, (I) 

f du dw\ dT dT , , NX 

w = 0, l—+e—\ +Re— =0, and — + Bi(T - T C (X)) = 0, on z = 1, (2) 
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where T C (X) = ln(l + <^X)/ln(l + <fr) is the conductive temperature profile in the cylindrical 
geometry and Bi is the Biot number. The Reynolds number Re = 7(^2 — T\)h2 / 'v2p^L is related 
to the Marangoni number by Ma = RePr where Pr = v/k is the Prandtl number of the fluid. 

Far from the lateral boundaries, the basic state (uq,0,wo,To) is determined as an expansion 
into the small parameter e, as performed by 0]. Up to first order in e, we obtain: 

uo(X, z) = Gv<S>{X)u {W, z) + £Gr 2 $p03ui(VF, z) + 0{e 2 ) 
w (X,z) = <D(e 2 ) (3) 
T {X, z) = T C (X) + eRa$(X)2T (W, z) + 0{e 2 ) 

with the Grashof number: 

r Lv2 

and the Rayleigh number Ra = GrPr. The parameter W — Re/Gr, which behaves like h~ 2 
is the inverse of the Bond number, it measures the relative importance of thermocapillary and 
thermogravity effects. The functions Uo(W, z) and To(W,z) are polynomials in z similar to those 
that describe the velocity and temperature profiles in rectangular geometry |2] and ui(W, z) is given 
in [7]. The main difference with the rectangular geometry is the non uniformity of the flow along 
the radial direction which appears through the slowly varying function §{X). In the forthcoming 
stability analysis this will prevent decomposition of the disturbances in Fourier modes in the radial 
direction. To bypass this difficulty we shall perform a local stability analysis and introduce the 
local Grashof number Grx = Gr$(X) as the governing parameter and similarly Rax = Ra$(X). 
In the following we shall assume that Grx = 0(1). Thus, the second term in the expansion of 
uq and To is of the order s$(X) — T(X) which is a small quantity only near the outer cylinder 
where r(l) e. Near the inner cylinder, T(X) takes values as large as T(0) = 0.475 and cannot 
be neglected. It will be shown in the next section that the perturbations equations depend on X 
exclusively through Grx and 



3 Linear stability analysis 

The basic state is perturbed by the superposition of three-dimensional disturbances v, for the 
velocity and T, for the temperature such as 

{v(r, 9, z, t),T(r, 9, z, t)} — expi(m9 + Lot) |w(r, z), v(r, z),w(r, z),T(r, z)| (4) 

where m is the azimuthal wave number and uj the frequency. The disturbances are assumed to 
evolve in the horizontal plane on a smaller scale than the basic flow. Thus, the fast variable 
x = e~ 1 X is introduced and the azimuthal wave number m is replaced by the local wave number 
13 = T(X)m. If the stability analysis is performed while T is kept constant, the disturbances can 
be sought as periodic functions of x such as for example u(r, z) — u(z) expiax. The linearized 
evolution equations for the perturbations reduce to a set of coupled differential equations 

Grx [((ict — T)uq + TGrx(ia — 3r)?2i) u 
+ (u' +TGt x u' 1 )w] 
Gix{ia + T)(u + TGr x ui)v 

iaGrx^o + rGrxUi)u> 



= -iap +{L-iu- T2)u - 2i[3Tv (5) 
= -if3p +(L-iu- T2)v + 2i(3Tu, (6) 

= -^- + (L-iu)w + T, (7) 
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Ra x [ia(u a + TGixu^T + (1 - 2TR& x T )u 

+ R& x fow] = (L-iPruj)T (8) 
with L — d2/dz2 — (a.2 + (32) + iaT. The incompressibility condition reads 

dw 

(ia + T)u + i(3v + — = 0, (9) 

The boundary conditions are 

u = v = w = T = 0, at z = 0, (10) 
io = 0, d z T + BiT = Q, d z u + iaWT = 0, d z v + i0WT = Q. at z = l. (11) 

When r = 0, the permutation of the cold and hot sides has no consequence on the instability 
pattern. Indeed, the transformation (Gr,a,u —> —Gr, —a, —u) leaves the above system invariant. 
This is no longer true when r ^ 0, in which case the sign of T needs also to be changed jointly for 
the system to be invariant. For fixed values of the curvature parameter T, the differential system is 
solved by the same method as in and critical values of the Grashof number, frequency and wave 
vector are found as functions of the curvature parameter T. 

The results presented on Fig.[21are obtained for values of the physical and geometrical parameters 
(Pr = 10, Bi = X,4>= 16, W = 6) corresponding to experiments described in UJ. On this figure, 
we have used the signed curvature sgn(T2 — Ti)T(X) and we can then observe the smooth behavior 
of all the quantities when crossing the zero value corresponding to the rectangular case. Results 
obtained for W = 2.4 [7j, which are not reported here, confirm the behavior of the critical quantities 
shown in Fig. 2. 



4 Discussion 

From Fig. 2c showing the critical value of the local Grashof number as function of the local curvature 
we deduced Fig. 2d showing the critical Grashof number Gr c versus the radial position, which is 
more convenient for comparison with experimental results. Particularly interesting is the decrease 
of the critical Grashof number when r — > 0. This means that hydrothcrmal waves will first appear 
near the center of the cylindrical cell, i.e., in region of higher curvature. An increase in the curvature 
also leads to a variation of the frequency (Fig. 2b) and over all to a significant variation of the wave 
vector orientation (Fig. 2a). 

When Ti > T\, the azimuthal wave number j3 vanishes for a value of the curvature T* = 0.24, 
which means that in cells with enough curvature, hydrothermal waves will be not only localized 
near the center, but they will propagate in the radial direction at onset. This fact was observed 
experimentally in 6] where shadowgraphic pictures of the instability pattern clearly show that at 
onset the spatial structure is made of concentric circles called targets, which are localized near the 
cold center of the cell. Slightly above onset the pattern spreads toward the outer side of the cell 
and deforms progressively to give rise to spirals that fill more and more space as the supercritical 
region extends (Fig. 2d), in agreement with experimental observations (Fig. 6,7 in [^]). 

When the center of the cell is the hot side (Tj > T2), the situation is quite different and there 
is experimental evidence that the instability pattern has always a non vanishing azimuthal wave 
number in agreement with our theoretical results. For a small height of fluid the pattern is made of 
spirals localized near the center of the cell jZj . For a higher height of fluid, an hexagonal structure 
is shown at onset 0. 
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5 Conclusion 



The flow which appears in a thin fluid layer confined between two differentially heated cylinders 
was considered. Far from the side walls, the basic state is modeled by a slowly varying flow along 
the radial direction. A local stability analysis is performed leading to the critical values of the 
Grashof number, wave numbers and frequency as functions of the local curvature. As a result, the 
instability is predicted to appear first near the inner cylinder and the values of the corresponding 
wave numbers are quite different whether the inner cylinder is the hot or the cold side. This is in 
agreement with experiments where spiraling waves have been observed near the hot inner side 
while pulsating targets appears when the inner cylinder is the cold side 
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Version francaise abregee 



1 Introduction 

Lorsque la surface libre d'une couche de fluide de faible epaisseur est soumise a un gradient ther- 
mique horizontal il apparait un ecoulement du a l'effet Marangoni qui resulte de la variation de 
la tension superficielle avec la temperature. Cet ecoulement qui est stationnaire pour de faibles 
valeurs du gradient thermique, se destabilise au-dela d'un seuil a partir duqucl il apparait des 
ondes propagatives hydrothermales selon la terminologie de Smith et Davis 0. 

Les analyses de stabilite [2] de ce type d'ecoulement etaient jusqu'a present surtout adaptees 
aux realisations experimentales en cellules rectangulaires |3J- Peu d'etudes theoriques traitent de la 
geometrie cylindrique 0] alors que de nombreuses experiences sont realisees dans cette configuration 
0- Pour y remedier, nous analysons l'effet de la courbure sur les proprietes critiques des 
ondes hydrothermales. 

Un schema du systeme considere est represents en figure f . Une couche de fluide d'epaisseur h est 
confinee lateralement par deux cylindres de rayons R\ et R2 = R\+L maintenus aux temperatures 
T\ et T2 respectivement. Le fluide est incompressible, de viscosite v, diffusivite thermique k, 
densite p, coefficient de dilatation thermique a. La tension superficielle a de l'interface varie avec 
la temperature T du fluide et da/dT = —7, avec 7 constant et positif. 

2 Etat de base 

Deux parametres sans dimension sont necessaires pour caracteriser completement la geometrie. 
D'une part, le rapport d'aspect e = h/L, et d'autre part le maximum de la courbure locale <fi = 
LjR\. Effcctuant le changement de variable radiale r = Ri + LX, la courbure locale L/r, est 
representee par la fonction $(X) = 4>/(l + <fiX). 

La vitesse et la temperature du fluide sont regies par les equations de Navier-Stokes dans 
l'approximation de Boussinesq et par Pequation de la chaleur. Un choix d'echelles approprie per- 
met d'ecrire ces equations pour des grandeurs sans dimension faisant ainsi apparaitre le nombre 
de Reynolds Re = j(T 2 - T 1 )h2/v2p Q L et le nombre de Grashof Gr = ag(T 2 - Tx)hi/Lv2 dont 
on notera le rapport W = Re/Gr. Apres multiplication par le nombre de Prandtl Pr = v/k du 
fluide on obtient respectivement le nombre de Marangoni Ma = RePr et le nombre de Rayleigh 
Ra = GrPr. Les conditions aux limites qui s'expriment par les relations £Q| et © representent 
une paroi inferieure rigide et conductrice et une surface superieure libre et plane ou les echanges 
thcrmiques sont decrits par l'intermediaire d'un nombre de Biot Bi. 

Loin des parois cylindriques qui continent le fluide lateralement, la vitessse et la temperature 
dans Petat de base (uo,0,wo,To) sont determinees par un developpement en puissance de s (pj). 
La variation lente de l'ecoulement de base dans la direction radiale est reportee dans la definition 
du nombre de Grashof local Grx = Gr§(X), permettant ainsi d'effectuer une analyse de stabilite 
locale ou Grx represente le parametre de controle et ou intervient le parametre geometrique r(X) = 
e$(X) mesurant la courbure locale. 
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3 Analyse de stabilite lineaire 



On superpose a l'etat de base des perturbations tri-dimensionnelles v pour la vitesse et T pour 
la temperature representees par des modes (@J de nombre d'onde azimutal m et de frequence uj. 
Nous considerons des perturbations qui evoluent dans le plan horizontal sur une cchcllc de longueur 
plus courte que ccllc de l'ccoulcmcnt de base, ce qui conduit a introduire la variable x = e X 
et a remplacer le nombre d'onde azimutal m par (3 = T(X)m. Pour une valeur fixee de T, les 
perturbations admettent comme solutions des fonctions periodiques en x, de nombre d'onde radial 
a. Apres linearisation, les equations d'evolution pour les perturbations se reduisent au systeme 
d'equations differentielles ©-©j avec la condition d'incompressibilitc © et les conditions aux 
limites associees ()10|l and (|1 1|) . 

En l'absence de courbure, T = 0, la permutation des bords chaud et froid n'a pas de consequence 
majeure puisque le systeme ©-© est invariant dans la transformation (Gr, a,u — > — Gr, —a, — u). 
II n'en va plus de meme lorsque r / 0, auquel cas il faut inverser conjointement F pour con- 
server l'invariance du systeme. Apres resolution numerique du syteme ©-© pour les valeurs des 
parametres (Pr — 10, Bi = l,<f> = 16, W — 6), on presente en figure 2 les valeurs critiques du 
nombre de Grashof, de la frequence et des nombres d'onde en fonction de la courbure. 

4 Discussion 

La variation de la valeur critique du nombre de Grashof en fonction de la position radiale est 
representee sur la figure 2d. On remarque une decroissance de Gr c lorsque r — * 0, revelant ainsi que 
Pinstabilite apparait en premier lieu pres du cylindre interieur ou la courbure est plus importante. 
Les variations de la frequence et des composantes radiale et azimutale du vecteur d'onde sont 
representees sur les Figs. 2a et 2b. Lorsque T\ < T2, le nombre d'onde azimutal s'annule pour une 
valeur de la courbure F* = 0.24 alors qu'il reste fini lorsque T\ > T2- Ces resultats sont en accord 
avec les observations experimentales effectuees par ombroscopie [HJEl. 

5 Conclusion 

Nous avons considere l'ecoulement qui apparait dans une mince couche de fluide confmee entre 
deux cylindres maintenus a des temperatures differentes. Loin des parois laterales, l'etat de base 
est modelise par un ecoulement lentement variable suivant la direction radiale. Une analyse locale 
de stabilite est effectuee conduisant aux valeurs critiques du nombre de Grashof, de la frequence 
et du vecteur d'onde en fonction de la courbure locale. Ces resultats montrent que l'instabilite 
apparait en premier pres du cylindre interieur avec un vecteur d'onde purement radial si ce cylindre 
est le bord froid et dans le cas contraire ce sont des ondes spirales qui apparaissent. Ces predictions 
theoriques sont en tres bon accord avec les resultats experimentaux exposes en [H], 0. 
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Figure 1: Section of the cylindrical geometry and associated notations. 
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Figure 2: Critical values of (a) the wavenumbcrs (a, f3), (b) the frequency u> and (c) the local 
Grashof number Grx as functions of the signed curvature, (d): Grashof number Gr versus signed 
radial position. 
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